TROLL model currently compute leaf lifespan with Reich’s allometry (Reich et al. 1991). But we have shown that Reich’s allometry is underestimating leaf lifespan for low LMA species. Moreover simulations estimated unrealistically low aboveground biomass for low LMA species. We assumed Reich’s allometry underestimation of leaf lifespan for low LMA species being the source of unrealistically low aboveground biomass inside TROLL simulations. We decided to find a better allometry with Wright et al. (2004) GLOPNET dataset.
We tested different models starting from complet model Mcomp: \[ {LL_s}_j \sim \mathcal{logN}(log({\mu_s}_j),\,\sigma)\,, ~~s=1,...,S_{=4}~, ~~j=1,...,n_s\]
\[{\mu_s}_j = {\beta_0}*e^{{\beta_1}_s*{LMA_s}_j^{{\beta_3}_s} - {\beta_2}_s*{Nmass_s}_j^{{\beta_4}_s}}\] \[{\beta_i}_s \sim \mathcal{N}({\beta_i},\,\sigma_i)\,^I\] \[(\beta_i, \sigma, \sigma_i) \sim \mathcal{\Gamma}(0.001,\,0.001)\,^{2I+1}\]
Figure 1: Leaf mass per area (LMA), leaf nitrogen content (Nmass) and leaflifespan (LL). Leaf mass per area (LMA in \(g.m^{-2}\)), leaf nitrogen content (Nmass, in \(mg.g^-1\)) and leaf lifespan (LL in \(months\)) are taken in GLOPNET dataset from Wright et al. (2004).
\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA + {\beta_2}_s*N,\sigma)\,\] Maximum likekihood of 10.0779021 and \(R^2\) of 12.653
\[ LL \sim \mathcal{logN}(\beta_0 + LMA^{{\beta_3}_s} + N^{{\beta_4}_s},\sigma)\,\] Maximum likekihood of -2.0074706 and \(R^2\) of 13.525
\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA^{{\beta_3}_s} + {\beta_2}_s*N^{{\beta_4}_s},\sigma)\,\] Maximum likekihood of 13.1931018 and \(R^2\) of 15.319
\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA + N,\sigma)\,\] Maximum likekihood of 1.0074785 and \(R^2\) of 12.877
\[ LL \sim \mathcal{logN}(\beta_0 + LMA^{{\beta_3}_s} + N^{{\beta_4}_s},\sigma)\,\] Maximum likekihood of -0.4205913 and \(R^2\) of 13.808
\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA^{{\beta_3}_s} + N^{{\beta_4}_s},\sigma)\,\] Maximum likekihood of 5.7451315 and \(R^2\) of 12.78
\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA,\sigma)\,\] Maximum likekihood of 4.2132211 and \(R^2\) of 15.128
\[ LL \sim \mathcal{logN}(\beta_0 + LMA^{{\beta_3}_s},\sigma)\,\] Maximum likekihood of -6.9667156 and \(R^2\) of 14.27
\[ LL \sim \mathcal{logN}(\beta_0 + {\beta_1}_s*LMA^{{\beta_3}_s},\sigma)\,\] Maximum likekihood of 8.787795 and \(R^2\) of 15.358
| M1 | M2 | M3 | M4 | M5 | M6 | M7 | M8 | M9 | |
|---|---|---|---|---|---|---|---|---|---|
| ML | 10.07800 | -2.00700 | 13.193 | 1.00700 | -0.42100 | 5.74500 | 4.21300 | -6.96700 | 8.78800 |
| R2 | 12.65282 | 13.52506 | 15.319 | 12.87726 | 13.80837 | 12.78011 | 15.12824 | 14.27014 | 15.35766 |
Figure 3: Model predictions.
Reich, P.B., Uhl, C., Walters, M.B. & Ellsworth, D.S. (1991). Leaf lifespan as a determinant of leaf structure and function among 23 amazonian tree species. Oecologia, 86, 16–24.
Wright, I.J., Reich, P.B., Westoby, M., Ackerly, D.D., Baruch, Z., Bongers, F., Cavender-Bares, J., Chapin, T., Cornelissen, J.H.C., Diemer, M. & Others. (2004). The worldwide leaf economics spectrum. Nature, 428, 821–827.